Resonant spin-changing collisions in spinor Fermi gases 
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Spin-changing collisions in trapped Fermi gases may acquire a resonant character due to the 
compensation of quadratic Zeeman effect and trap energy. These resonances are absent in spinor 
condensates and pseudo-spin- 1/2 Fermi gases, being a characteristic feature of high-spin Fermi 
gases that allows spinor physics at large magnetic fields. We analyze these resonances in detail for 
the case of lattice spinor fermions, showing that they permit to selectively target a spin-changing 
channel while suppressing all others. These resonances allow for the controlled creation of non-trivial 
quantum superpositions of many-particle states with entangled spin and trap degrees of freedom, 
which remarkably are magnetic-field insensitive. Finally, we show that the intersite tunneling may 
lead to a quantum phase transition described by an effective quantum fsing model. 



Cold spinor gases have recently attracted a large inter- 
est. A spinor gas is formed by atoms in two or more in- 
ternal states simultaneously confined by optical traps [l[ . 
Particularly intense efforts have been devoted to spinor 
Bose-Einstein condensates (BECs), which present a rich 
variety of ground state phases with different topologies 
0, 0, 0, SB], such as ferromagnetic, polar, uni- and bi- 
axial nematic, and more. Spinor Fermi gases are also at- 
tracting currently a growing attention. Collective modes 
in high-spin fermions have been shown to present rich 
features [7|. Spin-1 Fermions allow for color superfluid- 
ity and baryon formation, linking spinor fermion physics 
to QCD [8|. Spin-3/2 fermions also present fascinat- 
ing properties such as quintet Cooper pairing and Alice 
strings [9(. 

The dynamics of spinor BECs has been also actively in- 
vestigated, in particular the coherent oscillations between 
spinor components [lo| . In addition, spinor gases consti- 
tute a novel tool for the analysis of out-of-equilibrium 
systems, and the generation of topological defects after 
a rapid quenching across a phase transition [ll[ . On the 
contrary the spinor dynamics of high-spin fermions has 
been, to the best of our knowledge, not studied in detail. 

An external magnetic field typically constitutes a ma- 
jor handicap for the spinor dynamics. Although, due to 
spin conservation during a collision, the linear Zeeman 
effect (LZE) is irrelevant for the spinor dynamics, the 
quadratic Zeeman effect (QZE) may become (even for 
relatively low fields) sufficient to suppress spin change 
collisions, and hence any spinor dynamics However, 
under appropriate conditions an external magnetic field 
may even stimulate spin dynamics in spinor BECs, due 
to the compensation of QZE and mean-field shifts [12j . 

In the following, we show that spinor fermions allow for 
controllable magnetically-tuned resonances in the spin- 
changing dynamics, by compensating QZE and trap en- 
ergy. These resonances are typically absent in spinor 
BECs, either because the mean-field energy greatly over- 
whelms the trap energy in the Thomas-Fermi regime, or 
because in the weakly-interacting regime the trap levels 
are irrelevant since the bosons mainly occupy the ground 
state of the trap. In lattice arrangements, as those dis- 



cussed below, low temperature bosons occupy the low- 
est energy band, and hence the resonances play no role. 
They are also absent in pseudo-spin-1/2 fermions, since 
due to spin conservation the QZE energy is conserved. 
On the contrary, these resonances constitute a charac- 
teristic feature of high-spin trapped Fermi gases. 

In this Letter we illustrate the role of these resonances 
for the case of fermions with spin / = (2s+l)/2 (s > 1) in 
an optical lattice, where each site acts as an independent 
anharmonic trap, with up to three relevant levels (filling 
factor per spin component n < 3). However, we stress 
that the compensation between QZE and trap energy 
should also play an important role in more general opti- 
cal traps, leading to trap-dependent resonant spinor dy- 
namics, which will be the subject of a forthcoming work. 
Here, we study these resonances in detail for the lattice 
case, including interaction-induced shifts and resonance 
splittings, analyzing how for / > 5/2 these resonances 
may be employed to select particular spin-changing colli- 
sions while suppressing all others. Additionally, we show 
that an adiabatic sweep through the resonances allows 
for the controlled creation of quantum superpositions of 
states with entangled spin and trap degrees of freedom. 
These states are magnetic-field insensitive, and hence, 
contrary to the usual case, these systems allow for the 
study of spinor physics at large magnetic fields. Finally, 
we discuss the effect of the inter-site tunneling, showing 
that under proper conditions it may be described by a 
quantum Ising Hamiltonian, hence leading to a quantum 
phase transition between different multiparticle states as 
a function of the tunneling rate. 

The independent sites are described by the Hamilto- 
nian H = H + Hb + Hi. The trap energy is given by 



H 



n=0 



E(n) 



n=-f 



(1) 



where E(n) = hu>(n + (3n 2 ), and a n>m is the annihilation 
operator of fermions of spin m in the trap level n. The 
effective trap frequency (uj) and anharmonicity (/3) are 
obtained from the calculated first three on-site energies. 
The state of the system is described by a Fock state of 
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the form: \N 0i - f , ■ • • , N 2 , f ) ee ULo UL=-M, m Y 
where we keep the order (from left to right) from (n 
0, m = -f) to (n = 2, m = /) 0. 

The effects of a magnetic field £? are given by 
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where the LZE is given by p = g/isB, with g the Landc 
factor and /is the Bohr magneton, and the QZE by q — 
fJ,%B 2 1 'AhuJhs, where u>hs is the hyperfine splitting. 

Finally, the interatomic interactions (which we con- 
sider as being dominantly of short-range character) are 
provided by a Hamiltonian of the form |2| 
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where n = {711,712,713,714}, 771 = {mi, rrc.2, TO3, 7714}, 
C,i = / d 3 r<p n4 (r)*<p n3 (r^*<fi n2 (r)(f> ni (r), and C/ lR = 
J2F,M9F{f,mi,rn 2 \FM)(FM\f,rn3,mi). In these ex- 
pressions, gp = AmY?ap /m, with ap the scattering length 
for the collisional channel with total spin F (which due 
to symmetry must be an even number), M = —F, . . . , F, 
(/, mi, wi2 \FM) are the Clebsch-Gordan coefficients, and 
4> n (r) the 7i-th trap eigenfunction. The dynamics is 
calculated by integrating the corresponding many-body 
Schrodinger equation, using Runge-Kutta techniques. 

Although the resonances discussed below should be ob- 
servable for general initial conditions, the analysis of the 
effects implied is simplified by considering an initial mix- 
ture of 771 = ±1/2 (at the end of this Letter we discuss 
how this mixture may be achieved in on-going experi- 
ments). In addition, this initial state may be employed, 
as shown below, to controllably create complex superpo- 
sitions of spin-level entangled states. For a filling factor 
(per component) n = no ± 5n, with no = Q, 1, . . . , all lev- 
els up to Tio — 1 are filled in all sites, whereas the level no 
is occupied with a probability <5ti. This imprecise filling 
becomes eventually important, and is taken into account 
in our simulations. For n — 2, state |1) in Fig. [1] is the 
initial state for the dynamics. 

A collision leads to a QZE shift AEqze = q{m\ ± 
n\ — mf). As discussed above, this shift pre- 
vents spin-change collisions, as long as \AEq Z e \ is larger 
than the typical interaction energy of spin-change colli- 
sions. Level-changing collisions are characterized by an 
energy shift AE T rap = E (714) + E(n 3 ) - E(n 2 ) - E(ni). 
An interesting situation occurs when spin-changing col- 
lisions are at the same time level- changing. In that case, 
AEqze and AEtrap may compensate each other, lead- 
ing to resonances in the spin-changing collisions. Spin 
conservation precludes these resonances in usual spin-1/2 
Fermi systems, which thus do not show level-changing 
collisions if \AEtrap\ is larger than the spin-changing 
interaction energy. In the following, and for simplic- 
ity of our discussion, we concentrate in situations where 
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the harmonic trap energy can be assumed as conserved 
(m + 71 2 =713+ 714), and only the anharmonic change 
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is relevant. In general, 



AE. 

however, a sufficiently large B may lead to resonances 
which involve a violation of the previous condition, lead- 
ing to a wealth of resonances and complex multiparticle 
superpositions, similar to those discussed below. 
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Figure 1: States involved in the resonant dynamics at a ~ 
— 1/2 for / = 3/2. The lines indicate from bottom to top the 
three trap levels. The spin states are indicated by filled and 
hollow squares (m = ±1/2) and circles (m = ±3/2). 

The resonances arc altered by interparticle interac- 
tions. We illustrate this by considering the initial / = 3/2 
state |1) (see FigfT]). Note that multi-particle states |2) 
to |5) posses equal QZE plus anharmonic energy, Eqze + 
Eanh — 5q + 4:/3hu> (an important point as discussed be- 
low), whereas for |1) Eqze + Eanh = q + 2(3huj. Hence a 
resonance would be expected at a = — —1/2. How- 
ever, a detailed analysis shows that close to resonance 
the states indicated in Fig. [T] are isolated from the rest 
of all other possible (seven) states that can be reached 
via Hi from |1), leading to an effective 5-state problem. 
Interestingly, one may show that \D±) = (|3) + \2})/^/2 7 
and \D 2 ) = (|5) ± \i))/y/2 decouple from |1), playing a 
similar role as dark-states in quantum optics. The dy- 
namics is then fully dominated by the coupling of |1) 
and two "bright" states |±) = |[(|3) - |2)) ± (|5) - |4))]. 
As a consequence, the resonance splits into two peaks, 
corresponding to the energy shifts between |1) and |±): 



a R = -1/2 ± [ 5o £o + 92 (6 ± SO] /4/3^ (4) 

where £0 = C1100 — C2200 ± Clin/2, £1 = — 2Coooo ~ 
3C 2 2oo ± Cnoo ± Cim/2, and 6£ = 2C 22 oo- Note that, 
due to the EPR-like nature of the spin-change collisions, 
the resonances are accompanied by transitions into quan- 
tum superpositions of multi-particle states with entan- 
gled spin and trap level degrees of freedom, which we 
explore below. Split resonances for 1 < h < 2 (Fig. [5]) 
are hence a direct consequence of the coherent forma- 
tion of the above mentioned many-body entangled states. 
The split resonance peak may be employed to selectively 
detect |±) (which may become important to probe the 
linear superpositions discussed below), since initial |±) 
or |— ) states lead to different shifts of the resonance into 
|1), observable by monitoring the 771 = 1/2 population. 
The picture gets more complicated if 2 < fi < 3 (Fig. [2]), 
since the resonances are shifted depending whether there 
is none, one or two particles in the third trap level. Mul- 
tiple resonance peaks appear (or a resonance broadening 
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if the individual peaks are not resolved), but in this case 
due to the imprecise filling, as well as to the above men- 
tioned splitting. Since the spinor dynamics depends on 
n, one observes plateaux in d 2 N 3 / 2 /dn 2 (where N 3 / 2 is 
the sum of the time-averaged populations inm = ±3/2), 
with jumps when n becomes an integer (resembling the 
de Haas- van Alphcn susceptibility plateaux). 




Figure 2: (color online) / = 3/2. Time-averaged popula- 
tion (Nz/o) in m = ±3/2 as a function of a, for go, 2 = 
0.08, 0.10 [3, and n = 2.0 (thin, red), and 2.5 (thick). 
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Figure 3: (color online) / = 5/2. Time averaged popula- 
tion -/V3/2 in m = ±3/2 (thick) and N 5 / 2 in m = ±5/2 
(thin, red) as a function of a, for n = 2.0, and 30,2,4 = 
0.08,0.10,0.05 0. 

Whereas for / = 3/2 only one spin-changing collisional 
channel is possible, / > 5/2 allows for more channels. 
Fig. [3] shows the case of spin-5/2. Note that with an 
initial state |1) only two spin-changing collisions are pos- 
sible: from (ni 2 — l;™i.2 — =tl/2) either into (ri3 4 = 
0,2;m 3A = ±3/2) or (n 3A = 0,2;m 3A = ±5/2), with 
respective resonances at ~ —1/2 and cur ~ —1/6. In- 
terestingly, if the two resonances are well resolved, as it is 
the case in Fig. [3j it may be possible to select a particular 
spin-changing channel, preventing all others, providing a 
novel tool for the manipulation of spinor gases. As for 
the / = 3/2, the resonances couple the initial states with 
non trivial superpositions of spin-level entangled states. 

As previously mentioned, quantum superpositions of 
spin-level entangled states play an important role in the 
dynamics at the spin-changing resonances. These states 
may be created in a controllable way by sweeping adi- 
abatically the external magnetic field across the reso- 
nances. As an example, we consider the spin-3/2 case. If 



(3 < 0, the state |1) is the ground-state for \a\ > \an\. In 
that case, if g 2 > (32 < 0) an adiabatic sweep towards 
I a I < I or I creates |±) (|— )), with a fidelity larger than 
95% for sweeps with da/dt = 2.5 x lCT 4 ^. If the gas 
is prepared initially at a = or — 0.05, this linear ramp 
towards a = olr . ± 0.05 demands a change in 4000 trap 
periods, which for typical values of the effective on-site 
u) (few kHz) represents sweep times of less than 0.5s. 
Note that for (3 < 0, the created states are the ground- 
states for I a I < I or I, except for g 2 < 0,<?o > — 3^2 and 
92 > 0, .go > 52, where Z?i j2 are the (degenerate) ground 
states. However, as mentioned above, these states are 
dark states and hence they are not coupled to |1) during 
the sweep. Note that remarkably, all states belonging to 
the superpositions |±) possess the same QZE. This leads 
to two important consequences. On one hand, the sweep 
across the resonance is produced at rather large B (~ 1 
G, see below), but the result of the sweep crucially de- 
pends on the spinor physics (sign of g 2 ). Hence spinor 
physics becomes relevant even for large B. Second, once 
the states are produced an abrupt increase of B leaves 
the system unchanged, and hence the superpositions cre- 
ated may be robust at very large B, even opening the 
possibility of employing Feshbach resonances to variate 
92 or <7o 5 without polarizing the system due to the QZE. 
Similar states may be created by adiabatic sweeps for 
/ > 3/2. E.g., for / = 5/2 a sweep would create a state 
cos{4>) (1 2) - 1 3) ) / V2+ sinf» (j 4) - ] 5))/V2, where tan</> = 
(e - Ve 2 ±« 2 )M e = 15(2 52 - gi)(C Q220 - Gx>oo)/42, 
and n = 2(9g 4 - 2g 2 - 7go)C 1W 2/7. 

The states discussed above are created in isolated lat- 
tice sites. However, if the lattice is relaxed, the tunneling 
for the third band may become relevant, especially if the 
splitting 4g2Co220 between the states |±) becomes very 
small. In that case, different types of lattice Hamiltoni- 
ans may be created. Let us consider the particular case 
of spin-3/2, with \g 2 \ <C |go|- In that case, the states 
|±) form a quasidegenerate pseudo-spin- 1/2 doublet. If 
the tunneling t for the third band satisfies t \go\, then 
second order processes lead to an effective spin Hamil- 
tonian. Note that the condition \g 2 \ -C \go\ is neces- 
sary to avoid mixing other possible states. Employing 
l+X+l = 1/2+4*, I — >< — 1 = 1/2-3*, l+X-l = 
\ — )(±| = S x — iS y , the effective lattice Hamiltonian may 
be reduced to a quantum Ising Hamiltonian [lh| 



Heff = - Jeff SI 
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where J eff = 4:g 2 C 0220l and W eff = 
4t 2 2C 2222 -3Co22o _ Typically W eff > (< 0) if 

ogoC()220 C2222— 0O0220 J r J 11 K ' 

> (< 0). Hence there is a critical tunneling, 



*c ~ 3 32.9oC 2 2 2 ^ 



-3Cq2 



,....»„ 1 -77- rp 1, such that for t < t c 

— \ 202222— 0C0220 / 

the single-site physics dominates, and |±) (|— )) is 
the ground-state for g 2 > (g 2 < 0). For t > t c 
the tunneling dominates, and for go > the system 
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enters into a ferromagnetic state with all sites in either 
|u) = (|3) - |2))/V2 or \d) ee (|5) - |4))/V5. If g Q < the 
coupling is antiferromagnetic, and hence a ID system 
enters into a staggered configuration with alternated \u) 
and \d) states. The latter case opens the possibility for 
the analysis of frustration in other lattice geometries. 

Finally, we comment on possible experimental real- 
izations. We stress that the resonances should play 
an important role in high-spin Fermi gases in general 
dipole traps. For the case discussed above, the reso- 
nances should be clearly observable in standard lattices, 
Vq sin 2 kx. However, these lattices may be inappropriate 
for the detailed study of the rich phenomena discussed, 
since very large Vb > 50E rec (where E rec — ti 2 K 2 /2m) 
is needed to isolate the three lowest levels at each site 
from other sites, and it may be difficult to load only 
the first three bands of an initially weak lattice, due to 
the vanishing gap AE34 between the third and fourth 
band. A better scenario is provided by trimerized lattices 
(three wells per elementary cell) [IH . In our calculations 
we have considered for simplicity a ID trimerized lattice 
V(x) = Vo(sin 2 (Kx/3) + § cos 2 (2kiz;/3) + |sin 2 (Kx), as- 
suming a strong confinement in the yz-plane provided e.g. 
by an additional 2D lattice. V(x) allows for a relatively 
large AE34 even for relatively weak lattices. E.g. for 
Vo = 0.5£ rec , AE 3i ~ 0.22£ rec (for, e.g., ia K, typically 
Erec/ks — 2/iK), whereas the tunneling for the lowest 
bands is still sufficiently large (allowing for thcrmaliza- 
tion in the weak lattice). In this way, it may be possi- 
ble to prepare controllably a Fermi gas in the first three 
bands, if the temperature T < A_E 34 and n < 3. The con- 
sidered initial m = ±1/2 mixture may be created from 
a polarized Fermi gas in a maximally stretched state, 
by first employing an adiabatic radio-frequency sweep to 
transfer into m = 1/2, and then a 7r/2 radio-frequency 
pulse to establish a coherent m — ±1/2 mixture. A suffi- 
ciently large B (away from the discussed resonances) iso- 
lates this mixture against spin-change collisions due to 
QZE. In the presence of polarized bosons (e.g. in KRb 
mixtures), this mixture would thcrmalize in the weak lat- 
tice leading to an incoherent m = ±1/2 mixture, which 
after a lattice ramp-up would lead to the isolated sites 
discussed in this Letter. In a last stage the bosons should 
be eliminated, and B brought to the resonant values. For 
40 K, q/h ~ 70(f) 2 Hz, hence for 2tt/k = 800nm and 
V = 3.3E rec (w/2n ~ 9.25kHz, (3 = -0.3), a = -1/2 
occurs for B ~ 4.45G. 

Summarizing, spin-changing collisions in Fermi gases 
become resonant if the QZE compensates the trap en- 
ergy. These resonances are absent in BECs and spin-1/2 
Fermi gases, being a characteristic feature of high-spin 
fermions. We have shown for the case of lattice fermions 
that the resonances are shifted and split by interactions, 
and may permit to target a single spin-changing channel 
while avoiding the rest. The resonances allow for the con- 
trolled creation of quantum superpositions of states with 



entangled spin and trap-level degrees of freedom. These 
superpositions are magnetically isolated, and hence these 
systems allow for the observation of spinor physics at 
large magnetic fields, in principle even at Feshbach reso- 
nances. Finally, tunneling may compete with the on-site 
spinor physics to lead to a quantum phase transition de- 
scribed by a quantum Ising model. The rich phenomenol- 
ogy deriving from these resonances is by no means ex- 
hausted by the discussion above. Interesting physics is 
expected in the spinor dynamics of high-spin fermions 
in dipole traps. In addition, different resonances and 
spin mixtures in lattice fermions may allow for a wealth 
of quantum superpositions of spin-level entangled states, 
and various types of effective Hamiltonians. These per- 
spectives will be the focus of forthcoming works. 
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